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We present a detailed study of the magnetic structure and spin waves in the Fe jarosite compound 
KFe3(S04)2(OH)6 for the most general Hamiltonian involving one- and two-spin interactions which 
are allowed by symmetry. We compare the calculated spin-wave spectrum with the recent neutron 
scattering data of Matan et al. for various model Hamiltonians which include, in addition to isotropic 
Heisenberg exchange interactions between nearest (Ji) and next-nearest (J2) neighbors, single ion 
anisotropy and Dzyaloshinskii-Moriya (DM) interactions. We concluded that DM interactions are 
the dominant anisotropic interaction, which not only fits all the splittings in the spin-wave spectrum 
but also reproduces the small canting of the spins out of the Kagome plane. A brief discussion of 
how representation theory restricts the allowed magnetic structure is also given. 
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I. INTRODUCTION 

In search for unusual magnetic ground states and spin 
dynamics, frustrated systems have been the main fo- 
cus of both theoretical and experimental investigations 
in recent years i^i^i^iSiSiLSiSiiS In frustrated magnetic 
systemSji^ the energies of the various spin interactions 
compete and therefore they can not be simultaneously 
minimized. This competition can be induced by the ge- 
ometry of the lattice on which the spins are arranged, in 
which case the phenomenon is called geometrical frustra- 
tion and it often leads to quite unusual low temperature 
spin structures and dynamics. The nearest neighbor (nn) 
Heisenberg antiferromagnet (AF) on the corner-sharing 
Kagome lattice is one of the most studied systems since it 
has all the ingredients such as low dimensionality, strong 
frustration, and low coordination number, required for a 
disordered "spin Hquid" ground statei '^i'^i^i^i'^i^ 

Figure ^ shows three possible Neel states of the 
Kagome lattice nn antiferromagnet. In the classical limit, 
all ground states satisfy the "120° structure", in which 
the angle between each nn pair of spins is 120°. Fig.^i-b 
shows q=0 type ordering with respectively positive and 
negative chirality. For positive (negative) chirality the 
direction in which the spins rotate as one traverses clock- 
wise the triangle of sites is clockwise (counterclockwise) . 
Fig. ^ shows the a/S x a/S spin structure which has tri- 
angles with both positive and negative chirality. From 
the spin configurations shown in Fig. ^ it is clear the the 
classical ground state has a continuous degeneracy due 
to the "weathervane" rotation of the spins'^'"'"'''^'^'* and 
therefore no long range magnetic order is expected even 
at zero temperature. However, in real systems, there are 
small perturbations such as next nearest neighbor (nnn) 
interactionsffi' anisotropics or defects. It has been also 
shown that thermal or quantum fluctuations could also 
lift some of the continuous degeneracy known as "order 
by disorder," yielding coplanar spin structures 

Despite the extensive theoretical studies that suggest 




FIG. 1: (Color online) Three possible Neel states of the 
Kagome lattice nn AF: (a) q=0 with positive chirality, (b) 
q=0 with negative chirality, and (c) \/3 x \/3 structure. The 
plus/minus signs indicate the chirality of spins on the elemen- 
tal triangles. The shaded (yellow) areas indicate the magnetic 
unit cells. Note that the Neel sates shown above have contin- 
uous degeneracy as the spins on a line (a-b) and on a hexagon 
(c) can be rotated out of the plane as shown by ellipses with- 
out changing the classical ground state energy. 



many possible fascinating frustrated ground states for 
the Kagome nn AF, few experimental realizations exist 
for the system. The initial experimental studies were fo- 
cused on the layered garnet, SrCra;Gai2-xOi9 -' How- 
ever the interpretation of the magnetic properties of 
this system were complicated by the presence of an ad- 
ditional triangular lattice interposed between Kagome 
layers and by the inherent configurational randomness 
associated with the random alloying. In recent years, 
it has been shown that the jarosite family of miner- 



als AM3(OH)6 (804)2 fi2'i^'i^'iS*iLiLi2i22i2i (where A is 
monovalent ion such as K+, and M is trivalent cation 
such as Fe'^"'', Cr^+, or V^+) form much better reahzation 
of the two dimensional (2D) Kagome lattice. As shown 
in Fig. 12 for M=Fe jarosite, the magnetic ion M^~^ is 
centered in a slightly distorted and tilted oxygen octahe- 
dra, and it forms a Kagome lattice in the afe-plane. The 
Kagome planes are widely separated by nonmagnetic 
and SO^^ ions. With an exception of the hydronium 
jarosite (H30)Fe3(OH)g(S04)2, all members of jarosite 
family are found to exhibit long range magnetic order 
(LRO) at finite temperatures with varying ground state 
spin configurations and exchange interaction strengths 
depending on the magnetic ion M. The strength of ex- 
change interaction is the greatest for the parent M=Fe'^+ 
{d^,L = 0,S = 5/2) jarosite with Curie- Weiss tempera- 
ture Qcw = —800 K and T/v « 65 K. Replacement of the 
Fe^+ centers by Cr^+ (d^.S = 1/2) also affords an anti- 
ferromanetically ordered material but with a significant 
reduced Qcw — —67 K and Tjv « 2 K. Interestingly, 
when the magnetic ion is replaced by V^^, the ground 
state is changed to ferromagnetic Kagome layers which 
are coupled antiferromanetically."'^^ As we will discuss de- 
tail below, such a ferromagnetic ordering on the Kagome 
lattice is allowed by the representation theory in contrast 
to a previous analysis 

Among the members of the jarosite family, 
KFe3(OH)6 (804)2 (FeJ) is probably the most stud- 
ied one. The magnetic ground state of FeJ was 
first investigated by Townsend et al. using neutron 
diffractioniH However, later it became clear that the 
proposed spin configuration in Ref. |3 was not quite 
correct Inami et aii4 presented a detailed neutron 
scattering study and determined that the FeJ has 2D 
Kagome planes with the q=0 spin-structure with posi- 
tive chirality as shown in Fig. They observed that 
the unit cell along the c-axis is doubled in the magnetic 
phase, suggesting the the Kagome planes are antiferro- 
magnetically coupled. Field-dependent magnetization 
measurements^- have suggested that each Kagome plane 
has a small ferromagnetic component due to the canting 
of the spins in the q=0 structure with positive chirality 
(this is the so called "umbrella" configuration). However 
to the best of our knowledge, there is no experimental 
data which indicates the direction of this spin-canting 
whose determination would allow one to deduce the sign 
of some of the anisotropic terms in the Hamiltonian, as 
we will discuss in detail later. 

The fact that FeJ exhibits LRO at finite tempera- 
ture indicates that there are interactions other than nn 
isotropic AF superexchange, such as nnn interactions, 
anisotropics, etc. Inami et ali^ii^ have pointed out that 
in FeJ, the oxygen-octahedra are significantly tilted (see 
Fig. 121) , thereby inducing a strong single- ion crystal field 
(CF) anisotropy. They used this CF term in the Hamil- 
tonian to explain the observed spin configuration and to 
calculate the spin-wave spectrum. However in their spin- 
wave calculations, they neglected the canting of the spins 
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FIG. 2: (Color online) Conventional unit cell of FeJ. For 
clarity only Fe, big spheres (blue) and the octahedral oxygens, 
small spheres (red) are shown. The tilting angle, 60 of the 
FeOe octahedra is also shown. 



out of the Kagome plane. As we will discuss in detail be- 
low, this neglect leads to qualitatively wrong results for 
the spin-wave energy gaps at the high symmetry points 
in the Brillouin zone. 

Elhajal et alM- have pointed out that CF terms in FeJ 
should be fairly small because they are second order in 
the spin-orbit coupling and Fe'^"'" has a spherical charge 
distribution (z. e. L = 0). Hence as an alternative source 
of anisotropy, they suggested that the Dzyaloshinskii- 
Moriya (DM) interaction of the formiSiSi Di^Si x Sj 
could stabilize the experimentally observed spin struc- 
ture. Indeed, the z-component of the DM vector (the 
definition of its components is given below) forces the 
spins to lie in the a6— plane and therefore effectively acts 
like an easy-plane anisotropy. The sign of can distin- 
guish between the q=0 states with negative and positive 
chirality. 8ince the V3 x state has triangles with 
both positive and negative chirality, the Dz term alone 
will not affect the energy of this state. Furthermore if we 
introduce a Dy component, we break the rotational sym- 
metry around the c-axis and create a small anisotropy 
with respect to in-plane orientations. The effect of Dy 
is also to cant the spins (so that they have a small out- 
of-plane component) into the observed "umbrella" spin 
configuration. Finally, since the DM interaction occurs 
at first-order in the spin-orbit coupling, it is expected to 
be larger than the CF terms of Inami et alM^ 

In order to identify the origin of the magnetic interac- 
tions that are responsible for the LRO in FeJ compounds, 
clearly we need more experimental data such as the ob- 
served spin-wave spectrum. Fortunately, thanks to the 
very recent progress in the synthesis of high-quality sin- 
gle crystals of FeJ compounds^iSiSi such spin-wave data 
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has been recently become availablei^ Here we present a 
detailed theory for the spin-wave spectrum in FeJ com- 
pound for generic nn and nnn superexchange interactions 
including the CF and DM terms discussed above. As we 
shall see below, the DM term along with the nn and nnn 
isotropic interactions can explain quite well not only the 
observed spin configuration but also the observed spin- 
wave spectrum. 

Briefly this paper is organized as follows. In the next 
section we first discuss the symmetry of the FeJ structure 
in the paramagnetic phase and then present a brief repre- 
sentation analysis of the allowed magnetic structures. In 
doing so, we discovered that in Wills's analysisfi^ some 
of the ferromagnetic wavefunctions were missing. Indeed 
such a structure consisting of ferromagnetic easy-plane 
Kagome planes has been observed for a NaV-jarositeJ^ 
In this section, we also argue that characterization by 
irreducible representations is more fundamental than by 
the chirality of the q=0 state. In Sec. Ill we discuss the 
generic magnetic Hamiltonian that we use in our calcu- 
lations. In this section, we first discuss the symmetry 
of the Hamiltonian matrix and then its representation in 
a coordinate system in which the local z-axis coincides 
with the local direction of the spin moments. Here we 
also treat the canting of the spin orientations. In Sec. 
IV, we derive analytic results for the spin-wave energy 
gaps at the F, X and F-points of the Brillouin zone. 
In Sec. V, we present spin-wave spectra from numeri- 
cal calculations and compare the results with the recent 
spin-wave data of Matan et alM- Here we consider sev- 
eral models with increasing complexity. We find that 
the DM terms combined with nn and nnn isotropic ex- 
change interactions give an excellent fit to the data. Our 
conclusions are summarized in Sec. VI. The details of 
the diagonalization of the spin Hamiltonian were given 
in Appendices. In Appendix A, we describe the trans- 
formation to the boson representations of the spins. In 
Appendix B, we discuss how to obtain the normal mode 
energies and eigenvectors. Appendix C describes the ac- 
tual calculations of the matrix elements and Appendix 
D gives the results for the dynamical structure factor in 
terms of the normal mode energies and eigenvectors. 

II. CRYSTAL AND MAGNETIC STRUCTURE 

A. Symmetry of the Paramagnetic Phase 

We first discuss the symmetry of the paramagnetic 
phase. In Fig. |21 we show the conventional unit cell 
which contains three formula units of FeJ. The primi- 
tive unit cell (which contains one formula unit of FeJ) is 
rhombohedral, with basis lattice vectors 

ai = {a/2)i + (\/3a/6)j + (c/3)fc , 

a2 = -{a/2)i + {V3a/6)j + {c/3)k , 

aa = -(V3a/3)j + (c/3)fc , (1) 
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FIG. 3: (Color online.) Kagome planes. The filled circles are 
x-y planes at 2; = 0, the filled squares are x-y planes a.t z = 
c/3, and the filled triangles are x-y planes at z = 2c/3. The 
positive z-axis is out of the plane. Representative sublattice 
numbers r are given. 



where c denotes the height of the conventional unit cell 
of Fig. [5] The space group of FeJ is R3m, which is #166 
in the International Tables of Crystallography.^^ In order 
to understand the magnetic properties of FeJ we show in 
Fig. Oonly the Fe S" = 5/2 spins. Their locations within 
the primitive unit cell, denoted t„, are 

XI = (0,0,0) , T2 = (a/2,0,0) , 

X3 = (a/4,aV3/4,0) , (2) 

where the components are given with respect to the or- 
thogonal axes of Fig. 13 Notice that each plane (e. g. the 
filled circles) forms a Kagome plane lattice. As one moves 
from one plane to the adjacent plane at more positive 
the Kagome lattice is translated by {a/2)i + {\/3a/6)j, or, 
equivalently by either {—a/2)i+{\/3a/6)j or — (a\/3/3)j. 

Apart from translations the generators of the space 
group may be taken to be I, spatial inversion about the 
center of a hexagon, r, a two-fold rotation about an axis 
parallel to the x-axis and passing through the center of 
the hexagon, and R, a three-fold rotation about an axis 
passing through the center of the hexagon perpendicular 
to the Kagome plane. Since the center of a hexagon in 
one plane lies just above a triangle in an adjacent plane, 
the three-fold axis can be taken at the center of a trian- 
gle. The above operations imply the existence of mirror 
planes perpendicular to the Kagome plane and bisecting 
the sides of the hexagons. 

The reciprocal lattice basis vectors b^, which satisfy 
SLi ■ hj = 2'KSij, where Sij is the Kronecker delta, are 

bi/(27r) = {l/a)i + {l/V3a)3 + {l/c)k 
b2/(27r) = -(l/a)i + {l/v^a)j + (l/c)k 
b3/(2^) - -(2/V3a)j + (l/c)fc . (3) 

Since interactions between adjacent Kagome planes are 
very small, we may approximately describe the structure 
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FIG. 4: (Color online.) Reciprocal lattices for the FeJ system 
(blue squares) and for the two-dimensional hexagonal lattice 
(red circles) . The first Brillouin zone for the two-dimensional 
hexagonal system is indicated by heavy solid lines in red and 
that (in the 2 = plane) for the FeJ system by heavy dashed 
lines in blue. The light solid lines in red mark the boundaries 
of neighboring hexagonal zones. The FeJ Brillouin zone is 
three times as large as the hexagonal zone. 



in terms of a two-dimensional hexagonal lattice with basis 
vectors, A^, given by 

Ai = ai — a2 = az , 

A2 = a2 - ag = ^{a/2)i + (V3a/2)j . (4) 

The associated reciprocal lattice of this hexagonal lattice 
has basis vectors 

Bi/(27r) = (l/a)z + (1/V3a)j , 

B2/(27r) - {2/Vla)j . (5) 

In contrast the three-dimensional reciprocal lattice vec- 
tors in the z — Q plane for FeJ are 

bi - b2 = (47r/a)i 

b2 - bs = -(27r/a)i + {2'K^/i/a)j . (6) 

As shown in Fig. 01 the vectors of Eq. (O define the 
hexagonal Brillouin zones and those of Eq. © define the 
boundaries of the FeJ Brillouin zone in the z = Q plane. 
The first Brillouin zone of the FeJ system in the hexag- 
onal plane has an area three times that of the two- 
dimensional hexagonal system. If there were absolutely 
no interactions between Kagome planes, the spin-wave 
dispersion relation in the region in the first Brillouin 
zone of the FeJ system but outside that of the two- 
dimensional hexagonal system could be mapped into the 
dispersion relation inside the first Brillouin zone of the 
two-dimensional system and in this paper we adopt this 
picture. For the two dimensional system, spin-wave en- 
ergies at points just outside the hexagonal Brillouin zone 
near point Y are exactly the same as at the correspond- 
ing point just inside the zone near Y. For the three di- 
mensional FeJ system this is only approximately true. 



The spectra of these two points will differ due to inter- 
planar interactions. Similarly, for the two-dimensional 
system the spin-wave energies at point Z are identical to 
those at T, whereas for the FeJ system they will differ 
due to interplanar interactions. We will show the spin- 
wave spectrum observed by Matan et ah-- along the line 
T — X — Y\ where Y' is equivalent to Y if interplanar 
interactions are neglected. 



B. Results of Representation Theory 

A continuous phase transition may be described by a 
Landau expansion in powers of S7-(q), the Fourier com- 
ponents of the spin order parameters which are defined 
as 

s.(q) = (iyiv;;:)^(s(R-Kr))e*'i-(^+^) , (?) 

R 

where S(R-|-r) is the spin at site r in the unit cell at R, 
( ) denotes a thermal average, and q is the wavevector of 
the ordering. At quadratic order the Landau expansion 
for the free energy, F, assumes the form 

F ^ ^ Cr,Q;r',Q'(q)'5'r,Q(q)S'r',Q'(-q) ; (8) 
r,r';Q,a' 

where a labels Cartesian components and S't-,q(— q) = 
'S'r,Q(q)*- The ordering that actually occurs corresponds 
to the eigenvector associated with that eigenvalue of the 
quadratic form of Eq. (jSJ which first becomes negative 
(unstable) as the temperature is lowered starting from 
the paramagnetic phase. In view of the invariance of F 
with respect to the symmetry operations of the crystal, 
one can say that the eigenvectors give rise to an irre- 
ducible representation (irrep) of the space group of the 
crystali^ This analysis is straightforward and has been 
given previouslyfi^ although the results we present below 
seem to contradict that reference. 

Because spin is a pseudovector, spatial inversion takes 
a spin into its spatially inverted location but does not 
reorient the spin. Since all spins are located at centers 
of inversion symmetry, we only need to consider irreps 
which are invariant under spatial inversion, and we de- 
note these by r+, where the label n assumes the val- 
ues 1,2,3 and the superscript indicates invariance under 
spatial inversion. For the one dimensional irreps (F^ and 
rj) the situation is quite simple: each eigenvector of F is 
also an eigenvector of each operation of the space group. 
Thus we show in Fig. [Slthe possible spin configurations 
which correspond to one dimensional irreps. The actual 
spin configuration associated with irrep FJ is a linear 
combination of the two configurations shown. 

For the two dimensional irrep F^ the situation is more 

complicated. Here one has p wavefunctions ip)j , for 
i — l,2,...p which transform as the first row of the 
matrix representation of F;^ and their p partners V']'^'^^ 
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FIG. 5: Spin configurations which transform according to one 
dimensional irreps. Sohd circles indicate spins pointing out of 
the page. Top: a configuration of for which the eigenvalue 
of r is +1 and R is +1. Bottom: two configurations for Fj 
for which the eigenvalue of r is —1 and R is +1. 





which transform as the second row of the matrix repre- 
sentation of rj, where for FeJ p = 3^ We show these 
wavefunctions in Fig. El To be specific, if O is a symme- 
try operation and A4{0) is a matrix representation of it 
according to , then we have 



(3,«) 



FIG. 6: Spin configurations which transform according to 
two dimensional irrep F^. In the left (right) column we give 
the wavefunctions which transform according to the first (sec- 
ond) row of the matrices of the irrep. Arrows represent spins 
of unit length. Solid (open) circles indicate spins pointing 
perpendicularly into (out of) the Kagome plane. For vl/i^'^-' 
the magnitudes of the inward spins are both equal to half 
that of the outward spin which is set to be \/2. For ^'i'''^' the 



(9^ nonzero magnitudes of the spins are both •y/3/2 



These wavefunctions of correspond to the choice 



M{r) 



-1 
1 



MiR) 



-1/2 -V3/2 
V3/2 -1/2 



(ID) 



One can check that the wavefunctions shown in Fig. |^ 
actually do satisfy Eq. I^. Of course, all the wavefunc- 
tions shown here are orthogonal to one another and span 
the original subspace of spin components. (That appears 
not to be the case with the results of Ref . [iSj. In addition, 
his states do not seem to allow the net moment for a sin- 
gle Kagome layer to lie in the Kagome plane. Since this 
does happens for a Na jarosite^ the statement that the 
structure of the chromate jarosite must be the umbrella 
structure is not justified.) 



C. Discussion of Results: Dipolar Interactions, 
Chirality 

We now discuss the implications of the above results. 
The above results hold for two separate cases: in the fer- 
romagnetic case, all Kagome planes have the same spin 
structure (apart from a translation), whereas in the an- 
tiferromagnetic case successive Kagome plane have their 
spins inverted. The ferromagnetic case corresponds to 



wavevector q = 0, whereas the antiferromagnetic case 
corresponds to q = (37r/c)fcii^ Note that irrep r+ has 
two wavefunctions. This means that any such structure 
generically consists of a linear combination of a ferro- 
magnetic state (in a layer) with spins perpendicular to 
the plane and some in-plane ordering as well. This is the 
so-called "umbrella" state. The implication is that if one 
has an easy axis ferromagnet with spins perpendicular to 
the Kagome plane, one unavoidably will have some in- 
plane order. Conversely, in the case of FeJ one has order 
of irrep in each plane which must be accompanied by 
a net out-of plane moment. 

Note that the irreps and F^ are Ising-like in that 
within these representations one can not uniformly rotate 
the spins. Thus the excitation spectrum must have a sig- 
nificant energy gap. In contrast irrep is x-y-\ike. One 
spin structure is some linear combination of the wave- 
functions of the left column of Fig. Alternatively, the 
free energy is unchanged if one instead takes the same 
linear combination of the wavefunctions of the right col- 
umn of Fig. El The actual wavefunction is therefore 
characterized by four constants, a, /3, 7 and 9: 



tp = cos 9 aVP'i) /SV^'^^ 71/^^ 
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(3,2) 



-7-0, 



(3,2) 



(11) 



which can not be fixed by symmetry. At quadratic level 
in the Landau expansion, the free energy does not depend 
on 9 and one would expect a Goldstone (gapless) mode. 
However, due to the discrete symmetry of the Kagome 
plane, the anisotropy energy, F^^ which at lowest order 



Ea = Xcos(66') 



(12) 



will lead to a gap. If iiT > 0, the system will favor 6 = ir/Q 
(or equivalent angles) whereas ii K < Q the system will 
favor angles equivalent to 6* = 0. Quadratic energies, such 
as the dipolar energy, will not lead to a gap. Accordingly, 
it is not surprising that we found that for the dipolar 
energy of the system was independent of 9. 

It is interesting to note that if one assumes isotropic 
nearest neighbor (nn) and next-nearest neighbor (nnn) 

interactions, then spin configurations '^^^^ and ^'i^'' have 
the same energy. However, in general irrep will have 
lower energy because it can accommodate a distortion 

(2) 

into the configuration V&f, . In principle, local single- 
ion anisotropjiii and/or DM interactionsSS can select a 
ground state from among these functions. 

In the literature, there are frequent references to "chi- 
rality." Configurations have positive chirality, if when 
one traverses a triangle of spins clockwise, the spins ro- 
tate clockwise by A(/) = 120°, whereas for negative chi- 
rality A(/) — —120°. In our opinion, the characterization 
by irrcps is more fundamental. One sees that the wave- 
function of irrep and one of those of correspond 
to positive chirality and two of the wavefunctions for 
correspond to negative chirality. However, the sign of the 
chirality, of itself, does not indicate the presence or not 
of a gap in the excitation spectrum. Nor does the sign of 
the chirality correlate in an obvious way with the nature 
of other wavefunctions which may be admixed. 



D. Actual Magnetic Structure 

The actual magnetic structurc"'^"' of FeJ is that shown 
in Fig. 13 associated with the wavevector q = Stt/c. This 
means that successive Kagome planes are translated (as 
they would be in the paramagnetic phase), but then in 
successive planes the spin orientations are reversed. This 
causes a doubling of the unit cell because after three 
translations, the sites are back to their original positions, 
but the spin orientations are reversed. Thus the mag- 
netic unit cell contains twice as many layers as in the 
paramagnetic phase. Although each layer individually 
has a small ferromagnetic moment perpendicular to the 
Kagome plane, the sign of this moment alternates in sign 
from one Kagome net to the next, so that the sample ex- 
hibits no overall net moment. However, as in La2Cu04)i^ 
one has a spin-fiop transitionSi when a magnetic field per- 
pendicular to the Kagome plane rearranges the planes so 
that their ferromagnetic moments are parallel. 




^ ^ Spins in the second layer 



FIG. 7: Spin structure of FeJ as determined in Refs. IT^Hl . 
and|2ll 





FIG. 8: (Color online) Spin structure of FeJ. The squares 
(blue) are sites in the Kagome plane sd, z — c/3, the triangles 
(yellow) are sites in the plane ai z — —c/3 and the circles 
(red) are sites in the plane at z = 0. The two symmetry 
related choices for spin orientations in the ordered phase are 
shown in the two panels. 



Here we discuss some aspects of broken symmetry. One 
could ask whether the direction of the net moment of a 
Kagome plane is fixed by interactions within the plane 
or is accidentally selected. First, we should note that if 
we had a single isolated Kagome plane, then the two ori- 
entations or the net moment perpendicular to the plane 
would have identical free energy. As we will see, the di- 
rection of the net moment is fixed by the sign of the DM 
interaction. If we had a single Kagome plane, such an in- 
teraction would not be allowed. So, the uniqueness of the 
direction of the moment perpendicular to the plane, must 
depend on the details of the three-dimensional structure 
and, because +z and —z are equivalent in the param- 
agnetic crystal, on the in-plane magnetic ordering. To 
see how this is possible consider Fig. |S1 The actual spin 
structure is shown in the two panels of Fig. [S] The selec- 
tion of one of these states is arbitrary and represents an 
example of broken symmetry. Let us see how the direc- 
tion of the net moment perpendicular to the plane is fixed 
relative to the in-plane ordering. Suppose we have the 
broken symmetry selection of the structure shown in the 
left panel. There one sees that the in-plane spin moments 
point away from sites in the z = —c/3 plane and towards 



sites in the z = c/3 plane. This indicates that positive 
and negative z are not equivalent due to the presence of 
long-range magnetic order in the three-dimensional crys- 
tal structure. This means that when the spin order is as 
in the left panel, the free energy will select the direction 
of the perpendicular moment. When the spin order is 
as in the right panel, the perpendicular moment will be 
reversed relative to what it was in the left panel. As we 
will see, when the in-plane order is selected, the direction 
of the out-of plane moment is determined by the sign of 
the j/-component of the DM vector. It is interesting to 
note that the direction of the perpendicular moment can 
not be related to chirality. The structures shown in Fig. 
|S1 have positive chirality whether seen by a viewer from 
in front of the page or behind the page. So associating 
a direction with chirality can not distinguish between di- 
rections into or out of the page. 




FIG. 9: Numbering of sites used to define nn and nnn inter- 
actions. The axes fixed in the crystal are shown. 



III. MODEL HAMILTONIANS 



through the center of the bond. So the interaction matrix 
in Cartesian coordinates must be of the form 



A. General Hamiltonian 

We neglect interactions between adjacent Kagome 
planes. Accordingly, we now discuss the Hamiltonian 
of a single Kagome plane. We first parameterize the nn 
interaction between spins #1 and #2 in Fig. |51 which we 
write as 



Q/3 



(13) 



where the Greek superscripts label Cartesian compo- 
nents. The effect of the mirror plane perpendicularly bi- 
secting the 1-2 bond is to change the sign of the y- and z- 
components of spin and simultaneously interchange spin 
labels. Thus taking the transpose of the matrix and 
changing the signs of the y- and z-components must leave 
the matrix invariant. So the interaction matrix must be 
of the form 



M, 



5,2 



t(2) 



-D 



(2) 



D 



(2) 



t(2) 

t(2) 
'Jyz 



-D 



(2) 



y 
7(2) 

Jyz 



(2) 



(16) 



where D*^^) is the nnn DM interaction and J^^^ is the 
symmetric nnn interaction. 

We also include a single-ion anisotropy energy 
site i of the form 



for 



e: 



{l/2)Y,[C']r[STS^ + S^Sf] . (17) 

Q/3 



>i/]Q/3 r 



The mirror x-plane through site ^3 implies that the 
single-ion anisotropy matrix for this site is 















= 





c 

yy 


Cyz 


(18) 







C^'yz 


. 





Ml,: 



D 



(1) 



(1) 



(1) 7(1) 



J^y <Jyz 



yy 
(1) 



'Jyz 
^ zz 



, (14) 



where the superscript indicates an nn interaction and we 
introduce symmetric anisotropic exchange tensor J and 
antisymmetric DM interactionS(22i24 characterized by a 
DM vector D. Note that the definition of Eq. (^21 im- 
plies that 



(15) 



independent of what the local symmetry may be, where 
tilde indicates a transposed matrix. 

We next characterize the nnn interaction between sites 
#5 and #2 in Fig. |51 This interaction has to be invari- 
ant under the two-fold rotation about the x-axis passing 



B. Transformation to Local Uncanted Axes 

A direct, but inefficient, procedure would be to use the 
symmetry of the crystal to express all the other nn and 
nnn interactions in terms of the matrices of Eqs. 
and 1)16(1 . respectively. Instead, we proceed as follows. 
We express the nn and nnn interactions in terms of local 
axes which facilitate a spin-wave expansion. These local 
axes are defined so that the local positive z axis coincides 
with the projection of the local spin moment onto the 
Kagome plane. These axes are illustrated in Fig. [TUl 
The rotation matrices to transform to these local axes 



7^(ti) = 



1 



-i 

2 " 2 
1 



8 



nr2) = 



1 ^ 

1 

-10 
1 
1 



(19) 



The transformed interaction matrices, denoted I, are 
such that 



Q/3 



and 



n. 



52 



Q/3 



5,2 -^5 '-'2 5 



(20) 



(21) 



where now the spin components are taken in the local 
axes of Fig. ^Iso that 



Here 



2^1.2 
^5,2 



^(Ti)Mi,27e(T2) = J(l) 

7^(rl)M5,27^(T2)=X(2) 



-E. 



rjy Cjyz 



An) 
V 

(n) 



(") 



-E. 



(n) 



(22) 



(23) 



where 



£'(") 








~ 2 " 2 ^ 24 " 


1 


£;(") 
y 


- _ /(») 

~ '-'zz 

= <^n ~^^n ^" ^^^^ 




^ Z 


- 3 („) 1 („) V3 („) 

4^a: ^"^J/J/ 2 






2 2 24 


3 




- _in(") _ ^ /(«) 

- 2 y 2 






- n(") _ i /(") 






- /(") + 7(«)a _ in( 

— ^ K-'xx ^ -'yy ) 2^ 






- 2 2^ + 4 I 


A„ 
6 



2 



(24) 



where 



~ ^X£C '^zz "^yy ' ~ "^aia; "^zz 'V^^^ 




FIG. 10: Local axes of various sites. The local z-axis (indi- 
cated by the full line with arrow) is along the projection of the 
direction of the local magnetization onto the Kagome plane. 
The local y-axis is perpendicularly out of the Kagome plane. 
The local s-axis is in the Kagome plane and is indicated by 
the dashed line with arrow. 



so that J„ is the isotropic part of the nth neighbor in- 
teraction and we will only keep the anisotropic nn inter- 
actions, so that henceforth Aj = Ai and rjj = rji. We 
define i^i"' with a sign so that the interactions appear 
ferromagnetic in the local frame. The matrix elements 



(n) 

Eyz lead to a canted spin structure exhibiting weak fer- 
romagnetism, as we will see in the next subsection. The 
single-ion anisotropy in the local frame can be written as 



C = 
















Cyz 





Cyz 


Czz 



(26) 



A big advantage of expressing the interactions in terms 
of local axes is that all single-ion anisotropics and all nn 
interaction matrices are identical (apart from transposing 
if site indices are interchanged) . Using the three-fold axis 
and inversion symmetry we find that 



2^1,2 — 2^2,3 — 2^3,1 — 2^2,1 — Iz,2 — 2^1,3 
= 24^3 = X5,4 = X3,5 = X3,4 = 14^5 



1-5,3 



(27) 



where the sites are numbered as in Fig. 1101 

Just as for nearest neighbors, all the next nearest 
neighbor interaction matrices are the same, providing we 
take the indices in the correct order. 



C. Canting 

Up to now we have taken the local z-axes to lie in 
the Kagome plane. However, the matrix elements E^ 
for n = 1 or n = 2 gives rise to a field which induces 
a uniform y component of spin. For the spin- wave cal- 
culation it is convenient to define "canted local axes" 
such that the canted z-axes lie along the direction of the 
canted spins. This direction is found by minimizing the 
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FIG. 11: The 'initial' and 'final' local y and z axes. Initially 
the z axis is in the Kagome plane. The sign of 6 is the same 
that of Hy, which is the negative of -B^V , i .e. the negative of 



energy. We write the ground state energy per site, Eq^ 
for Sz — S cos 9 and Sy = sin as 

Eg/S'^ = [Czz - 2EP - 2Ef^] cos^ 9 

+ [Cy,-24i)-242)]3in2(? 

+2[Cyz + 2£;W + 2£;g)] sin 6* COS 6* 
= -Acc cos^ 9 - Ass sin^ 9 + 2 Acs sin 9 cos 9 



--[Acc + Ass] - HzCos{2e) 
-Hysm{29) , 



where 



Hy = -Ac 



(28) 
(29) 



[Acc — Ass]/2 , 
The energy is minimized by setting 

cos(26') = Hz/H , sin(26l) = Hy/H , (30) 



where H = ^H^+Hy. If we neglect the effect of 

anisotropic nnn interactions on the canting angle, we 
have 

Hy = -[Cyz + 2E^^}] = -Cyz - VSDy + Jyz 



Hz = E^) - '-Cz 

+eP - 42) 

= '-J.-^Dz 
2 2 



Cy-z 



:Cy 



where Ca-p = Caa — Cpp. Here and below the super- 
script "1" for nn's is implied in the anisotropic interac- 
tions. 



as illustrated in Fig. ^] where c = cos 9 and s = sin 9. 
Then, if we include this transformation we finally arrive 
at the local anisotropy matrix as 














^yy 


Cyz 





Cyz 


C zz 



c = 



and the transformed local interaction matrices are 



(33) 



2^1,2 = 



J,y 



-d-Z -Ey ^ E^y] 

dy Ey^ -E^ 



(34) 



and similarly for X5,2, where 

(?CzZ + S^Cyy -f 2csCyz , 



C zz 
in) 



E 
C 

E. 

C 



yy 
■(») 
y 



Czz T' L- ^yy ^^^^yz ; 



c (^1111 2iCsC^^j 



Eyz - 

dz = 



Cyz COS(20) + ]^[Cyy " C z z] ^^{29) , 

M';) cos(26l) + i - M")] sin(26l) 
,v - 2 



y 



(35) 



An equivalent condition for equilibrium is that 

'^^yz - (36) 



Cyz + 2YE[:^=0. 



When the x-z elements are summed over all neighbors 
they give a vanishing effective field. Indeed, group theory 
guarantees that this is the case for general interactions. 



IV. CALCULATION OF THE SPIN WAVE 
SPECTRUM 

In Appendix ^ we write the Hamiltonian in terms 
of boson operators and in Appendix ^ we obtain 
the wavevector dependent spin-wave dynamical matrix 
(SWDM) M(q) whose eigenvalues are the spin-wave en- 
ergies at wavevector q and which is written as 



M(q) = 



A(q) 


-B(q) 


B(q)* 


-A(q)* 



(37) 



This means that we want to transform the initial local 
coordinates with moments in the Kagome plane into the 
canted local coordinates via 



~*loc,i 



;?loc,i 



-rloc.f 



y 



(32) 



A. General Form of the Spin- Wave Dynamical 
Matrix 

The matrices A(q) and B(q) are given in Eqs. HA7|I 
and ljA8|l . respectively. We write them as 



A(q) = A£ + ^A(")(q) 



(38) 
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and 



form 



B(q) = - Cyy)£ + J2 B^"' (q) , (39) 



where A = C'xx 



yy 



2Czz, £■ is the unit matrix, and 



for n = 1 or 2 A^"' (q) is 





-2idi"'l7l"2'(q) 


+2i4"']7["3'(q) 


r ??(") 

[—-Cjx — -C/j, 

+2idi"']7{"i(q) 




-24"']7^:3'(q) 


-2^di"']7l"3'(q) 


+24"V^"3'(q) 





,(40) 



and B(")(q) is 






I-eL-^+e';'] 


i-Ei"^ + e';']- 






X7("3'(q) 









X7l"2'(q) 




X7S(q) 









x7S(q) 


X7i"3'(q) 





,(41) 



where 7^"^-^ (q) is the normahzed form factor for the nth 
shell of neighbors. Here we we wiU only treat n ~ 1,2. 
For in-plane interactions we have 



7^3^ (q) 
7f.^(q) 

7^3' (q) 



co8{aqx/2) , 
cos[a(q:r + V3qy)/A] , 
cos[a((j^ - V3qy)/A] , 
cos{aqyVi/2) , 
cos[a(3qa; - V3qy)/4:] , 
cos[a(3qa; + \/3gy)/4] . 



(42) 



In Eqs. 1001) and (|^ one one uses Eq. to relate the 
overlined coefficients to the bare coefficients, which are 
defined in Eqs. ^ and if^ . 

To simplify the expressions we will now specialize to 
the case when the nnn interactions are isotropic and all 
interactions further than nnn are ignored. Then we write 
the matrix A(q) in the form 





{Ai^iai)y[lh 
{^J2-ia2)-y[f 


{Ai+iai)yg^+' 
(iJ2 -ia2) 7^3' 


{Ai+iai)^[l^ + 
(^2 + ^2)7^2' 


Ao + 2J2 


{Ai^iai)yg^ + 
(ij2 -^2)723' 


{A^-za,)'yg^ + 
_ (iJ2-ia2)7w 


{Ai+iai)'y^23 + 
(|J2 +iQ2) 723' 


Ao + 2 J2 



(43) 



where 7!™ denotes 7,m'i(q). The matrix B(q) is of the 



Bn 


-^1 112 \^/ 

-|^27f2'(q) 


B-i'y'iV (a) 

-^1 /IS 

-|J27l3'(q) 


-|J27f2'(q) 


Bo 


Bi7^s'(q) 
-|^27g'(q) 


Bi7l3'(q) 
-1^271? (q) 


Bi-/g\ci) 

-f^27f3'(q) 


Bo 



(44) 



In Appendix lUl we evaluate of the constants in the ma- 
trices for the general Hamiltonian introduced above. In 
the interest of simplicity we now limit consideration to 
the following generic Hamiltonian: 



JiSi ■ Sj -f Dij ■ Si X Sj 



+ ^ J2Sk • Si 



DY,{sf'r-EY,[{sfr-{sfn (45) 



where £ nn (e 
over nn's (nnn's). 



nnn) indicates that the sum is 
Aj = {0,Dy{i,j),Dz{i,j)) is the 
Dzyaloskinskii-Moriya vector for bond i — j as shown in 
Fig. Inland the single-ion anisotropy terms are those used 
by Nishiyama et alr-^ in their treatment of the spin- wave 
spectrum in jarosites. Here the prime-spin components 
refer to the local axis associated with the rotated oxy- 
gen octahedra shown in Fig. 1 (see Ref. for details). 
The Dy{i,j) and D^(i, j) are all expressible in terms of 
the parameters Dy and Dz as discussed in Sec. HI. The 
single-ion anisotropy constants D and E used in Ref. 
are related to our generic single-ion matrix C defined in 
the uncanted-local frame given in Eq. 126|l as 



Cxx — 

Cxy — 

Cyz — 

Czz 



/ 

CxZ ^ Cyx C zx 

i:>cos(e'o)^ - £;sin(6lo)^ 

C,y = {D + E) sin(^o) cos(0o) 

i:'sin(6'o)^ - EcoaiOaf , 



(46) 



where 6q is the rotation angle of the octahedra around 
Fe-ion (see Fig. El and is 6*0 = 20° i& 

In what follows we will often have recourse to two sim- 
ple models. In the first of these, which we call the "DM" 
model, we neglect the single ion anisotropy and the ex- 
change anisotropy, so that the only nonzero parameters 
are Ji, J2^ Dy, and Dz- In the second model, which we 
call the "CF" model, all the anisotropy is incorporated 
by the single-ion crystal field, so that the only nonzero 
parameters are Ji, J2, and Cap- In either case, we as- 
sume Ji to be the dominant interaction. Accordingly we 
give here the expressions for the matrix elements of the 
SWDM which pertain to these two cases: 



Ao = 2Ji + A- 2V3Dz 



— [3DI + 2CI] 
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Bo 
Bi 
ai 
a2 



2 2 



1 



ByJ2 



6Ji 

2^/lDyD, DyJ^ 



Jl 



where A = Cxx + Cy 



3Ji 



2a 



(47) 



B. Analytic Results for Zero Wavevector 



When a is nonzero, we split the degeneracy between 
the two heretofore degenerate frequencies. This spht- 
ting is hnear in a. In contrast, when a = 0, the spht- 
ting between this two-fold degenerate mode and the other 
mode win be found to be of order \/ Ji5, where 5 involves 
anisotropic exchange or DM interactions. 

Note that for C = 0, we have stability for < D^c, 
where 



J 



V3 / 



Q \ J 



0[{Dy/jf 



(53) 



which is similar to Elhajal et atfS, who give D^c/J ~ 
Q.22{D,Jjf. For D = 0, we have stability (in the large 
J limit) if both Cxx — Czz and Cyy — Czz are positive, so 
that the z-axis is really the easiest axis. 



At zero wavevector the SWDM, M, whose eigenvalues 
give {^1 S) = Cj, is 



^0 


Ai - ^a. 


Ai + .a 


-Bo 


-Bl 


-Bl 


Ai + .o 


Ao 


Ai - ^a 


~~Bi 


-Bo 


-Bl 


Ai - iQ 


Ai +iQ 


Ao 


-Bl 


-Bl 


-Bo 


Bo 


Bl 


Bl 


-Ao 


-Ai - ia 


-Ai + ia 


Bl 


Bq 


Bl 


-Ai + ic 


-Ao 


-Ai - 


Bl 


Bl 


Bo 


-Ai - ia 


-Ai + ia 


-Ao . 



where a — ai + a2 and we include the effects of J2 by 
replacing Ji in Eq. 147|l everywhere by Ji + J2. Take |1) 
to have components (1, 1, 1, 0, 0, 0)/\/3 and |2) to have 
components (0, 0, 0, 1, 1, l)/\/3. Then 

M|l) = (Ao + 2Ai)|l) + (Bo + 2Bi)|2) 

M|2) = (-Bo-2Bi)|l) + (-Ao-2Ai)|2) . (49) 

From this we find the spin-wave energy loq to be 

((I'd)" = (Ao + 2Aif - {Bo + 2Bif . (50) 

Now take |3) to have components 

(2, — 1, — 1, 0, 0, 0)/\/6, |4) to have components 

(0, 0, 0, 2, — 1, — 1)/V6, |5) to have components 

(0, 1, — 1, 0, 0, 0)/\/2, and |6) to have components 
(0,0,0, 0,1, -1)/V2. In this subspace 



M = 



a b i\/3a 

—b —a i^/ia 

-i^/ia a b 

—i^/3a —b —a 



(51) 



where a = Ao — Ai, b — Bo — Bi, and a = ai + a2. 
Thereby we find that the other two spin-wave energies, 
, are given by 



uj± — \/ a?' — b"^ ± \fZc 



(52) 



Thus for zero wavevector we find the frequencies a) for the 
DM and CF models to be those given in Table n (under 
the heading "With Canting"). 



C. Results for Other Wavevectors 

1. The X Point 

The X point is at 

Qx = [47r/(3a)] , Qy = Qz ^ 



and the form factors are 
1 

"2 ' 
and 



11,2 



^(1) 

'1,3 



7i,2 — 



(2) 1 
72,3 = -1 



7i,3 — 



Thus the SWDM is of the form 



ail 


ai2 


-ai2 


-bii 


-bi2 


bi2 


ai2 


0.11 


-ai2 


-bi2 


-fell 


bi2 


— 0,12 


-ai2 


an 


bi2 


fel2 


-fell 


bu 


bi2 


-bx2 


— Oil 


* 

-0-12 


ai2 


bi2 


611 


-bi2 


—012 


—ail 


ffll2 


-bv2 


— bi2 


bii 


ai2 


ai2 


— Oil 



where, to the same accuracy as before 

an = ^0 + 2 J2 , 

= 2{Ji + J2)-2V3Bz + A 
+2Dl/J^+AClz/{iJi) , 



ai2 = -[Ai - iai]/2 + J2/2 - ia2 , 



-Clzl{l2Ji) - 



3i DyJ2 



2 Jl 
, 2{DyDz/Ji)]/6 



(54) 



(55) 



(56) 



(57) 



(58) 



n2 



''12 



(59) 



12 



TABLE L Spin wave energies uj = uj/S for zero wavevector for the DM model (only Ji, J2, Dy, and nonzero) and for the 
CF model (only Ji, J2 and Cap nonzero). The results neglecting canting are discussed in subsection IVE. Here J = Ji + J2. 
In this and succeeding tables Ca-p denotes Caa — Cp/s. 



With Canting 




DM Model 


CF Model 


Wo 


VT2\Dy\+0{Dyj) 
ADyD./J + 0{D'^/.f) 


[12 JC,^. + 4C,_.C,_, + 4Cl] + 0(C5/V J3/2) 
[6JC„_, + 4a_.C,_, + 7C,2,] + 0(C5/VJ»/2) 
2Cy. + O(CVJ) 


2 

Cj+ — Cj- 


Ncglecting Canting 




DM Model 


CF Model 






2V3Dy 


\yiJ(Jx — z ~\~ ^Cx — zCy — z] ^ 
[QJCy — z ~\~ ^Cx — zCy — z] ^ 




2 

U1+ — OJ- 



TABLE 11: Spin wave energies uj = uj/S at the X point for the DM and CF models. Here J2 is considered to be of the same 
order as D or C. The results neglecting canting are discussed in Sec. IVE. 



With Canting 




DM Model 


CF Model 


[Do 

2 

— Ul- 


[6Ji(3J2 - VSD,) + 3Dl - 

^[Ji + J2] - + 0{D^/J) 

{2DyD,-3s/3DyJ2) 1 (^(-£)3yj2-| 


[6Ji(3J2 + Cy-,) + 4a_,Cy_, + 7C^, 

+i2j2a-.]i/2 + o(c5/Vj3/2) 

^[Ji + J2] + %2C^-. + C^-.) + o(cVJ) 
Cj,. + O(CVJ) 


Neglecting Canting 




DM Model 


CF Model 


2 

[I)+ — Cj- 


[6Ji(3J2 - VSD,) - + 18Dl] ^'"^ 

^[Ji + J2] - ^i>. + ©(dVJ^) 

^Dy 


[6Ji(3J2 + Cj,-.) + 4a_,Cj,-. + 1272^-.]'/' 
^[Ji + J2] + ^(2C:.-. + Cj,-.) + o(cVj) 




foil 



— — Cy y — 26*^2/ (3 Jl) 



(60) 



and 



612 = -Bi/2 - (3/2) J2 



-J2 



V3 



4Ji 



/(12Ji) 



(61) 



where we have neglected symmetric anisotropic exchange 
and only kept terms that seem to be most relevant. 

Now we diagonalize the submatrices a, a*, and b. The 
eigenvectors of these 3x3 matrices are 



1 

7! 



A 
1 

-A2 



(62) 



where A'^ = 1, so that Ai = 1, A2 = (—1 + iA/3)/2, and 
A3 = (— 1 — i\/3)/2. One can show that the vector given 



in Eq. (|62|l is an eigenvector of a with the associated 
eigenvalue 



a(A) = ail + A^ai2 + Aa^j 



(63) 



and is simultaneously an eigenvector of b with eigenvalue 
6(A) = 611 + \Hi2 + A612 . (64) 
The three spin- wave energies are given by 



a(A) + a(A*) 



a(A) -a(A*) 



1/2 



(65) 



so that 

cZ^o = [(an + 2a'i2)' - (foil + 2foi2)Y^' 



[(an 



(foii-foi2f]'/'±\/3a'i'2 .(66) 



We thus find the spin-wave energies to be those given in 
Table ini 
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2. The Y Point 

The Y point is at qy — 27r/ (\/3a), Qx = Qz = 0, so that 

7l"i = 0, 7^"3^ = 0, 7i'i = 1, and -f^^l = -1. Therefore 
the SWDM at the Y point is of the form 



where 



On 


ai2 





-foil 


-bi2 





a*i2 


an 





-fol2 


-foil 











an 








-foil 


bii 


bi2 





-an 







bi2 


fell 





-ai2 


-an 











foil 








-an 



(67) 



The perturbed (i. e. for a'h 0) secular equation is 

;^L°)f]p^-(4"'f] + (a'/,)4 

,// N2r~,2 I 2 I J \2 1 2 , i2 







-2[a'l^f[Cj' + al^-{a',^f-b\,+bl^] 
Thereby we find the approximate solutions 

(a'1'2)' 



UJ± 



~(0) 



± 



{uff+al,-{a',2f-bl,+bl. 



(75) 



(76) 



In Table IIIII wc give analytical results correct only to 

leading order in the anisotropics, in which case uj± = 

~(o) 

10+ . 



ai2 



ail 



= ^ + 2 J2 , 

= 2(Ji + J2) - 2%/3I?, 



A 



+2DllJi+WlJ{iJi) , 



[^1 
(^1 



iai\ 



^2)/2 - 



£12 

y 



6J1 



(68) 



- J2/2 + ia2 , 

2 2Ji 

i^/i[Cy, + 2{DyD,/Ji) - 2i?„ J2/ Ji]/3 



il2 



il2 



foil — — Cj:a; — Cyy — 2 C^^ / (3 Jl) 



(69) 



(70) 



and 



fol2 



Si + 



A) 



2 



3Ai 



6Ji 



(71) 



One root for {ujY (coming from rows and columns ^^3 
and #6) is 



fo?i 



(72) 



and the remaining eigenvalues come from the projection 
of the SWDM, M, into the remaining subspace: 



V. FIT TO THE EXPERIMENTAL SPECTRUM 

Recently Matan et al^ have measured spin- wave spec- 
trum in FeJ using inelastic neutron scattering and here 
we discuss the comparison of our calculations with their 
results. The experimental spectrumJS. shown in Fig. ^] 
clearly indicates three spin-bands. One of these bands 
is quite dispersionless and it is therefore identified to 
be the "zero-energy" mode, which is lifted to about 8 
meV due to anisotropic magnetic interactions which we 
wish to identify. In addition to this "lifted zero-energy 
mode", the spectrum shows two additional gaps at 1.7 
meV (an in-plane mode) and around 7 meV (an out- 
of-plane mode). The spin- wave spectrum extends up to 
20 meV at the X-point. It is of interest to develop a 
model Hamiltonian which reproduces the observed spec- 
trum. In particular, we will see that the observed energy 
splittings at the various high symmetry points provide 
sensitive constraints on the interaction parameters. 

In order to identify magnetic interactions in FeJ, in 
this section, we will present numerical calculations of the 
spin-wave spectrum for a large number of models and 
compare the results with the experimental data shown 
in Fig. 1121 For a quantitative comparison we define an 
error-factor R (i.e. the goodness-of-fit) as follows: 



R 



100 

^data 



E 



[qi] - Wcai(qi)l/^oxp('7j) (77) 
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-foi2 


ai2 


an 
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-an 


-a*i2 


. fol2 


foil 


-ai2 


-an _ 



(73) 



As usual, in our analytic work, we work to second order 
in the anisotropics, in the absence of which ai2 = 012 + 
ia'h is real. In that case we have the unperturbed mode 
energies 



(an±a'i2)'-(foii±foi2)' 



where A'data is the number of experimental data points, 
and Woxp(9i) and Wcai(9i) are respectively the experimen- 
tal and calculated spin wave energies at wave- vector qi. 
One should recognize that R describes the percentage er- 
ror averaged over the entire spectrum. It is important 
to reproduce the splittings at the high symmetry points 
even though this only weakly affects the value of R. 

We start by comparing the experimental spectrum to a 
simple model where we consider only the isotropic nearest 
neighbor interaction Ji and a non-zero in DM inter- 
action vector (i.e. = Dy = 0). The best fit from this 
(74) simple two-parameter model is shown in Fig. 1131 The DM 
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TABLE III; Spin wave energies u = u)/ S a.t the Y point for the DM and CF models. The results neglecting canting are discussed 
in Sec. IVE. 



With Canting 




DM Model 


CF Model 




[l2Ji J2 - eVSJiD, + 4J| - 


[l2Jl J2 + QJlCy-z + 4Ca;-zCj;-z + 8J2Ci-z 




+18DI + 3Dl] + 0(D5/V^'^') 


+2J2Cj;-z + 4J| + 7C^,] + 0{C^/^/J'^/^) 




2Ji + f J2 - + 0(I>VJ') 


2Ji + 1 J2 + i(6a-z + 3C«-z) + O(CVJ) 


2 


Cj+ — (Do 


J2 + + OiD^J) 


J2 + a-z- iQ-z+C(CVJ) 


Neglecting Canting 




DM Model 


CF Model 


UJ- 


[l2Ji J2 - 6%/3JiD, + 4J| - + 18Dl] 


[l2Jl J2 + 6JlCy-z + 4Ci:_zCh-Z 

+8J2C:.-z + 2 JaCy-z + 4J|] 


2 


2 Ji + fJ2 - + 0(1)77) 


2Ji + f J2 + i(6C.-z + 3C«_z) + O(CVJ) 


Qi+ - Cjo 


J2 + + 0(p^l.]) 


J2 + C,_z- iCj,_z+0(CVJ) 
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FIG. 12: (Color online) Experimental spin-wave spec- 
trum in FeJ obtained from inelastic neutron scattering 



is discussed in the text in detail. 



interaction with only the Dz term lifts the zero-energy 
mode to the dispersionless energy S\/—6V^JDz, as ob- 
served experimentally. However this term alone can not 
give the experimentally observed small dispersion of this 
mode. Without Dy, we do not get any canting of spins 
either. The effect of the term is almost identical to 
an easy plane anisotropy. The "lifted zero-energy mode" 
is found to be degenerate with the other spin gap {i.e. 
the out-of-plane spin-gap) at the F point. However, we 
do not get the in-plane gap as experimentally observed 
at 1.7 meV. (Analytic results are given in Table D) 

The agreement between calculations and experiment 
can be improved significantly by turning on Dy in the 
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FIG. 13: (Color online) Spin-wave spectrum from a simple 
two-parameter J — _Dz model. In these and succeeding figures 
the values of the parameters of the fit are given in meV. 



DM-vector as shown in Fig. E| The main effect of Dy 
term is to create a small in-plane anisotropy and therefore 
it can be adjusted to give an in-plane gap of 1.7 meV, as 
experimentally observed. This model gives three gaps 
at r point, for which the analytic results are given in 
Table m From these analytic results, one sees that the in- 
plane gap is proportional to \Dy \ and the small splitting 
between the flat-mode and the out-of-plane gap is equal 
to lODyDz/ J , which is quite small. 

Although the experiment does not give clear results for 
the splittings of the high-energy modes at the X and Y 
points, we see from Tablellllthat the splitting, 5, at the X 
point is expected to be quite small (it is second order in 
the perturbations). Table Hill indicates that the splitting 
of the high frequency modes at the Y point is expected 
to be larger than at the X point. 

The Dy term also causes the spins to be canted with a 
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FIG. 14: (Color online) Spin-wave spectrum from a simple 
three-parameter Ji — Dy — model. 



FIG. 15; (Color online) Spin-wave spectrum from isotropic 
nn and nnn interactions Ji and J2. 



canting angle of about 2.5°. However despite a significant 
improvement with the addition of Dy term, we still do 
not get any dispersion to the flat mode in contradiction 
to the experimental finding. We also studied the single 
ion term D and E as given in Eq. H4()|l with only Ji 
and obtained similar results. With only DM or single 
ion anisotropy without further neighbor interactions, it 
was not possible to obtain the small dispersion of the flat 
mode as experimentally observed. 

In order to get the experimentally observed dispersion 
along with the correct average energy of the nearly flat 
mode, one needs to consider next-nearest neighbor inter- 
actions. In Fig. 1151 we show the best fit from a simple 
isotropic nn (Ji) and nnn (J2) interactions, which is in 
good agreement with the previous work.'^ We note that 
the nnn J2 interaction gives the zero-energy mode a sig- 
nificant dispersion. However in order to correctly repro- 
duce the flat-mode energy requires a J2 which gives too 
strong a dispersion. Since the Hamiltonian with isotropic 
Ji and J2 is rotationally invariant, we do not get any gap 
at the r point, in contrast to experimentally observed 
three gaps. 

From above discussion, it is clear that one needs both 
a nnn J2 interaction and either a DM or a single ion 
anisotropy. In Fig. ^| we show the best fit from a 
model which has the isotropic nn Ji , nnn J2 interactions 
and the DM vector {0, Dy, Dz). The fit to the data is 
very good. We reproduce not only the three gaps but 
also the small dispersion of the flat-mode. The spin gaps 
at the r, X, and Y points obtained numerically are in 
good agreement with the analytic results given in Tables 
m and mil The spin-canting angle is about 2.1°, in 
reasonable agreement with the experimentally estimated 
value of 1.5° M 

As an alternative to DM interactions, one can also try 
to explain the experimental data based on a model which 
considers only the single ion anisotropics D and E as 
given in Eq. (|46|l .^^ In Fig.lltib. we show the best fit from 



such a model which has an isotropic nn Ji and nnn J2 
and the {D, E) single ion anisotropy instead of the DM 
interactions. The R- factor for this fit is slightly larger 
than the R-factor obtained from DM interaction. The 
main differences are the larger splittings of the optical 
modes at the F and X points and the smaller canting 
angle for the CF model in comparison to the DM model. 
The canting angle from the CF model is only 0.8°, which 
is smaller than the experimental value^^ of 1.5°. 

The biggest problem with the CF model is the large 
gaps at the F and X-points as shown in Fig. Il()b . The 
experimental splitting of the modes at these points are 
much smaller. Numerical results give a splitting of about 
0.75 meV at the F point and 0.4 meV at the X-point. 
The splitting at the Y-point is quite small. These results 
are in good agreement with the analytical results given 
in Tables nil and Iml 

In conclusion, the observed spin-wave spectrum can be 
explained by a simple model which has only nearest and 
next-nearest isotropic interactions. The DM interaction 
gives the best fit to explain the observed gaps and the 
spin canting. Even though the single ion terms also fit 
the data well, the quality of the fit is not as good as the 
one from DM term at the F and X points. For the Fe ion, 
the single ion anisotropy is expected to be smaller than 
the DM term (because the single-ion term is second order 
in the spin-orbit coupling while the DM term appears in 
first order). For these reasons, we think that DM term is 
a natural source of anisotropy in FeJ compound, which 
not only stabilizes the observed experimental spin struc- 
ture but also gives the right spin- wave spectrum with the 
proper gaps and dispersion. 
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FIG. 16: (Color online) Spin-wave spectrum from isotropic nn 
and nnn interactions Ji and J2 and (a) DM- vector (0, Dy, D^) 
and (b) single ion anisotropy terms D and E. 



A. Effect of Spin-Canting on the Spin-wave 
Spectrum 

Here we give results when the transformation to the 
canted structure is omitted, i.e. 9 is forced to be zero. 
One purpose of giving these results is to show the impor- 
tance of making this transformation when calculating the 
gaps in the spectrum which depend on small perturba- 
tions. Also, we could check our calculations by comparing 
these results to those of Nishiyama et allS- who omitted 
this transformation. When one distorts a coUinear spin 
structure by the application of a small transverse field, 
the effects on the spectrum are second order in the dis- 
tortion angle. Here, however, the situation is different 
because the spins before distortion are noncoUinear. As 
a result of this noncoUinearity the matrix element d^,^'' 
in Eq. H23|l (which would be zero for a coUinear system) 

is of order J. Thus the matrix element c?^ in Eg . H35|l 
has a contribution, sd^"'' , induced by canting of the same 
order as that, C(ii"\ which existed before canting. This 
fact indicates that canting should not be ignored in the 
calculation of the spin- wave spectrum. 



To obtain results when canting is ignored we replace 
quantities with overbars in Eq. (|35|l by the corresponding 
quantities without overbars. 



1. Zero wavevector 

Here (where J denotes Ji + J2) 

Ao = 4£;^ + A = 2 J - + A , (78) 

Ai = -E, -Ey^ J 12 + ViD,l2 , (79) 
ai = 24(1) = -Dy , a2 = , (80) 

= CxX — Cyy , (81) 

Bi = Ex- Ey^ -3 J/2 + V3D^/2 . (82) 

Then, when the anisotropy in J„ is neglected, we find the 
frequencies listed in Tabled under the heading "Neglect- 
ing Canting." 

2. The X point 

The treatment including canting is now modified so 
that 



2( Ji + J2) - 2V3D, + A , 

ai2 = -[Ai - iai]/2 + J2/2 - ia2 

= -Ji/i + J2/2-V3D,/A-iDy/2, 

-B1I2 - 3J2/2 

3J1/4 - 3J2/2 - %/3Dj4 . (83) 



an 



hi2 



Then we find the spin-wave frequencies as listed in Table 

nn 



3. The Y point 



Now we have 



an 


= An(q 


= 0) = 2(Ji + J2) - 


ai2 


= ^12 (q 


= 0) = (Ji - J2)/2 


foil 




Cyy , 


hi2 


= -3(Ji 


- J2)/2 + ^/iD,/2 



(84) 

so that the frequencies are those given in Table ITTTl Un- 
less the only anisotropy is due to Dy, these results are 
the same as those given in Table lHll with canting, so that, 
except for this case, the effect of canting at the Y point 
is very small. 
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4- Numerical Results and Discussion 

The effect of canting can be deduced from the epxres- 
sions for the matrix elemets given in Appendix C. One 
sees that the rotation to the canted spin structure in- 
duces terms of the form D^/Ji into the SWDM matrix 

elements. Also the quantity ai = 2c?^ ^ is transformed 
(by the rotation to the canted system of coordinates) 
from the value Dy (neglecting the anisotropy in J„) to 
{Cyz + 2DzDy/ Ji)y/3- DyJ2/Ji. S upr pr Islugly perhaps, 
this rotation through an angle of order Dy/Ji also intro- 
duces a contribution to 012 of order J2Dy/Ji. A qualita- 
tive understanding of what is going on here is obtained 
by noting that the coefficient dz is associated with an 
interaction proportional to S^°'^'^ 3^°'^'^ , where the com- 
ponents are taken in the "final" i. e. canted coordinate 
system. This term directly affects the spin-wave energies. 
This term arises from rotating an interaction in the lo- 
cal uncanted coordinate system of the form dyS^°'''^S^°'^'\ 
a term which for coUinear systems is zero, but for this 
noncoUinear system is of order Ji, as given in Eq. (|24|l . 
Since these quantities control several of the splittings, it 
is clear that neglecting canting can have a major impact 
on these splittings, as shown in Fig. 1171 and the tables. 
From Fig. I16f a). it is clear that the quality of the fit is 
significantly degraded for the DM term, in particular at 
the r and X-points where big gaps are opened. 

The effect of neglecting canting is opposite for the 
CF term. The quality of the fit is improved because 
the large splitting at F (see Fig. llBr b)) is now zero 
(Fig. Elb)). For instance, for zero wavevector the split- 
ting a)+ — a)_ (given by 2\/3a) is linear in Dy when cant- 
ing is neglected, whereas when the canting was taken 
into account the splitting was found to be proportional to 
Cyz+2DzDy/ Ji. To calculate this splitting it is therefore 
not qualitatively correct to ignore the canting. In addi- 
tion, when canting is ignored and single-ion anisotropy 
is not present, neglecting canting leads to ljo = 0, even 
when D is nonzero. 

At the X point the effects of canting are similar. For 
the low energy mode the transformation to the canted 
structure introduces a term proportional to Dy which is 
only important if all the other anisotropics are very small 
or vanishing. However, as Table^indicates, the splitting 
between uj^ — uj- is proportional to Dy when canting 
is neglected, but is proportional to Cyz + 2DyDz/Ji — 
3y/3J2Dy/ Ji when canting is taken into account. 

At the Y point the effects of canting are less important. 
Indeed loq and are insensitive to canting and only in 
uj^ does one see that canting introduces a term of order 
D^ 

In conclusion, it is important that one should not ig- 
nore the spin-canting when spin- wave spectrum is calcu- 
lated. Even though the CF single ion term gives much 
better fit when the canting is ignored (see Figs. El and 
I17|l . it is not an appropriate approximation, as we have 
shown. When the CF single ion term is properly treated. 
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FIG. 17: (Color online) Spin-wave spectrum for the DM (a) 
and CF (b) models when the spin-canting is ignored. We used 
the same parameters as in Fig. 1161 to emphasize the effect of 
spin-canting on the spectrum. 



it actually gives a worse fit to the data and therefore it 
is more logical to conclude that the dominant source of 
the anisotropy in FcJ compound is the DM term. 



VI. CONCLUSIONS 

We have presented a detailed study of the mag- 
netic structure and spin-wave spectrum in FeJ for the 
most generic nn and nnn exchange interactions includ- 
ing Dzyaloshinskii-Moriya and the single ion anisotropy 
terms of Inami et aZ.^^ Our conclusions can be summa- 
rized as follows. 

• We first discuss the allowed magnetic structures 
from a representations analysis. The experimen- 
tally observed "umbrella" structure of q=0 with 
positive chirality arises from a one dimensional 
representation which exhibits three-fold symmetry. 
Our analysis and discussion of broken symmetry 
suggests that the characterization by irreducible 
representations is more fundamental than that by 
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chirality, even though chirahty has been widely in- 
voked in the htcrature. 

Representation analysis reveals the existence of an 
x-y-\ike two dimensional representation which has 
sixth order in-plane anisotropy. In this representa- 
tion one can have Kagome layers with a net in-plane 
moment which can be stacked either ferromagneti- 
cally or antiferromagnetically. The analysis of Ref. 
ITsI did not find wavefunctions of this type. How- 
ever, the observed magnetic ground state of a NaV- 
jarosite^* consists of such ferromagnetic Kagome 
planes which are coupled antiferromagnetically. 

We presented detailed spin-wave calculations for 
the most generic spin Hamiltonian allowed by sym- 
metry. We derived analytic expressions for the 
spin- wave energies at the T, X, and Y points in 
the Brillouin zone. The agreement between our an- 
alytic and numerical results tends to support the 
correctness of both approaches. Since, the spec- 
trum is invariant under a change of sign of Dy , the 
sign of this parameter remains undetermined by ex- 
periment. 

We obtained the spin-wave spectrum numerically 
throughout the Brillouin zone for a large number 
of models with increasing complexity and compared 
the results with the recent experimental data.^^ We 
obtained a nn exchange constant of Ji = 3.225 
meV, which is in good agreement with previous 
measurements We found that the observed small 
dispersion can be explained by a nnn isotropic in- 
teraction of J2 = 0.11. The observed gaps at the 
r, X, and Y points are best explained by consid- 
ering only the DM interaction with \Dy\ = 0.218 
meV and Dz = —0.195 meV together with nn and 
nnn interactions. The term is adjusted to lift 
the zero-energy mode to about 8 meV with a small 
dispersion due to J2 interaction. A fit of similar 
quality can also be obtained using the crystal field 
(CF) model of Ref. 1^ but the resulting anisotropy 
is about 10% of Ji, too large a value for the S state 
Fe'^"'" ion which is supposed to be quite isotropic. 
The calculated canting angle is 2.07° and 0.77° for 
DM and CF models, respectively, compared to the 
experimentally deduced value-^^ of 1.5°. 

Finally we discuss the effect of the canting on the 
spin- wave quantization axis which was neglected by 
Inami et al.^^ We found that neglecting canting is 
not a good approximation for calculating the spin- 
wave energy gaps at the high symmetry points in 
the Brillouin zone. When the CF is treated prop- 
erly, it gives larger gaps at the F point than those 
from DM interaction. The experiment does not in- 
dicate significant splitting at the F point and there- 
fore it suggests that the DM term works better to 
explain the data. 
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APPENDIX A: BOSON HAMILTONIAN 

To construct the boson Hamiltonian it is convenient to 
introduce spherical components {+^x + iy,~' — x~ iy) 
for interactions between site r in unit cell i located at Ri 
and site r' in unit cell j located at Rj . If these sites are 
nth nearest neighbors, then we write 

'-^TT'i}^i) ~ '^[-^xx ^'^-^xy ~f-^yx ~^-^yy] i-^^) 

and similarly for the single- ion anisotropy matrix C. Now 
we invoke the Holstein-Primakoff transformation in the 
form 

S^, ^ S- ala,r , S+ = V2Sa„ , = V2^a|(A2) 

Then apart from a constant the Hamiltonian in terms of 
Bose operators is 



2SY. 



-ST.' 



(A3) 



where the prime on the summation excludes the term 
with (i,r) = {j,T'). Now write 



so that 



(A5) 



where Eq is the ground state energy and 
niq)/S = ^A^.(q)at(q)a.(q) 
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-^^^BafiiqTaaiqjaf^i-q) , (A6) 



a/3 

where 



-xX(q)+irr^(-q) 



(A7) 



= [C,^ - Cyy]Sr,r' + 2j;;, i-q) , (A8) 



where 



T"f,(q) = ^T"f,(z,,)e--(-+^-^--') , (A9) 



F"./3 / ■ ■ 



where R.y = — Rj and the superscripts assume the 
values + , — , and z. Since each site is a center of inversion 
symmetry, the exponential factor in the sum over i in Eq. 
[9)1 can be replaced by a cosine: 

T"f,(q) = ^X,"f,(z,j)cos[q.(r + R.,-Tl]AlO) 



and I(q) is an even function of q. If G is a vector of the 
reciprocal lattice (such that G • R is a multiple of 27r), 
then 

Arr'{q + G) = A,,,(q)e^°-(^'-^) 

= ^,,,(q)cos[G. (x'-r)] 

= B^^/(q) cos[G • (t'-t)] . (All) 

Here we used the fact that for the Kagome lattice, the 
t's are half a lattice vector, so that e'*^ '''' ^'''^ = cos[G • 
(T'-r)]=±l. 

APPENDIX B: NORMAL MODES 

We write the normal mode operators as 

^l(q) = EK(q)«a(q) + <(q)««(-q)] , (bi) 



and these are determined by 

[■H,Xl{q)]^J'^\q)Xl{q) 



(B2) 



If there are p spins per unit cell, then we expect that p of 
the uj'--^^ will be nonnegative and p will be nonpositive. 
Equation HB2|) gives 



E 

Q,/3 



{q)Y,K{q)ai{q) + d^{q)a^{-q)] 

a 

c£:(q)^/?,a(q)a|3(q) + c^ici)Blp{q)a(,{-q) 
di^{q)Aaj3{-q)ap{-q) 

(B3) 



E 

a,/3 



+<(q)Sa,/3(-q)4(q 

where again u = u/S. 

1. Eigenvalue Problem 

This leads to the eigenvalue problem 



(q)c^:(q) = 5]K.Mq)c^(q)" 5/3,0 (-q)d^(q)] 
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(q)C(q) = Y.^B0,^{qYc^p{q)~ApA-^)d>;{m) 



which can be written as 

A(q) -B(-q) 
B(q)* -A(-q) 





■ c''(q) ■ 




■ c''(q) " 




.d^(q). 


.d^(q). 



(P5) 



where c'^(q) is a column vector with components c^'(q), 
C2 (q), . . ., c^(q), where p = 3 is the number of spins in 
the unit cell. 

We now note that because the Hamiltonian is Hermi- 
tian and taking account of inversion symmetry, we have 

Arr'iq) = A,v(q)* =^rr'(-q) . (B6) 

Similarly we write 

Brr'{c\) - B,,,(-q) =B,,,(-q) . (B7) 

Thus the eigenvalue problem is 



A(q) -B(q) 
B(q)* -A(q)* 



cf(q) 
d^(q) 



cD(^)(q) 



c'^(q) 
d^(q) 



(B8) 



We now show that the roots come in pairs with oppo- 
site signs. Equation l|B8|) is 



Ac'' - Bd'' 
B*c^ A*d^ 



(B9) 
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Take the complex conjugate of these equations, change 
the signs of both sides of each equations, and reorder the 
equations to get 

A[d'']* -B[c^]* = -tl.(^)[d'']* 



so that 



B*[d^]* - A*[c^]* = 



K]* , (BIO) 



where we used the fact that the a)'s are reaL This shows 
that from an eigenvector with eigenvalue tl;'^'^) we can 
construct a related eigenvector with eigenvalue — ti;^''-'. 
One can show that the eigenvalues at wavevector q are 
identical to those at wavevector — q, as expected in view 
of the inversion symmetry of the lattice. Thus we may 
write 



A(q) -B(q) 
B(q)* -A(q)' 

c(q) d(q)* 
d(q) c(q)* 



c(q) d(q)' 
d(q) c(q)* 



-cD(q) 



(Bll) 



where c(q) and d(q) are matrices whose columns are the 
column vectors c^(q) and d'^(q), respectively, and uj is 
a diagonal matrix with entries a)'^^-'(q). 



2. Dependence on Reciprocal Lattice Vector 

Recall the dependence of the matrices on the recipro- 
cal lattice vector as recorded in Eq. (|A11(I . Then the 
eigenvalue problem at wavevector q + G may be written 
in the form 

a;(^)(q + G)c^(q + G) = ^[A„,^(q)c^(q + G) 

-i?/3,„(-q)d^(q + G)]e'«-(^^-^°) 
c:-(^)(q + G)d'^{q + G) = ^[i3^^„(q)*c^(q + G) 

-A0,„(-q)d^(q + G)]e^°-(^''-^°) . (B12) 
Set (q + G) = cZi^^) (q) and 



c^(q + G) = c;^(q)e-'°-^'' 
d^(q+G) = d^(q) 



-iGT 



(B13) 



Then we recover the previous equations. So we con- 
clude that under addition of G, the eigenvalues remain 
invariant, but the wavefunctions change according to Eq. 



3. Orthogonality 

Taking the Hermitian conjugate of eigenvalue equa- 
tions, we get 

c''(q)tA(q)-d''(q)tB(q)t = (::,(^)(q)c^(q)t 
c''(q)^B(q) -d^(q)tA(q) ^ u''^'^^)d''{qy (BU) 



_c^(q)t ,~d^(q)t 
= Cj^^'Hq) \ c^(q)t ,-d^(q)t 



A(q) -B(q) 
B(q)* -A(q)* 



(B15) 



Thus from any right solution we have constructed the 
associated left solution. Now consider 



W^, ^ [c'^(q)t ,-d^(q)t] 



A(q) -B(q) 
B(q)* -A(q)* 



(B16) 



Using the fact that the left and right vectors are eigen- 
vectors, we have 



Wf,, = cZ'^'^Hq)* [c^(q)t ,-d^(q)t 



c-(q) 

d^q) 



(q) c'^(q)t ,-dA'(q)t 



d-{q) 



(B17) 



If /i = 1/ this shows that (I'('^)(q) is real. Also 



c'^(q)t ,-d^(q)t 



d-{q) 



0, fi^iy . (Bli 



(If we have degenerate eigenvalues, we can choose the 
wavefunctions to preserve this orthogonality.) 
A more general orthogonality is 



ct(q) -dt(q) 
-d(q) c(q) 



c(q) d(q)' 
d(q) c(q)' 



J , (B19) 



where X is the 2p x 2p unit matrix. 

4. Summary 

We have the transformation to normal modes as 

At(q) = ^[cj:(q)at(q)+<(q)a„(-q)] ,(B20) 

a 

and also 

A^(q) = J2K{qraUcD+di^iqral{-q)] (B21) 

a 

where fi is restricted so that a;^'^^ > 0. We require that 
Xfj^iq) obey Bose commutation relations: 



[X^{q),Xl{q)]^S^,, 



(B22) 



This gives 

5](cif(q)*<(q)-di!(q)*<(q)) = S,,. . (B23) 
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The inverse transformation is 

a„(q) = ^[cj:(q)X^(q)-<(q)*X4-q)t] (B24) 

We record here the dependence on reciprocal lattice 
vector: 

al(q + G) = ^[c^(q)*X^(q)t-d^(q)X^(-q)]e^«-^° , 
a^q + G) = ^[c^(q)X,,(q)-dj:(q)*X4-q)t]e-'^ff26) 



To use this we would assume that q is in the first Bril- 
louin zone and G is the vector needed to bring the actual 
wavevector back into the first Brillouin zone. 



APPENDIX C: EXPLICIT EVALUATION OF 
MATRIX ELEMENTS 

We now develop explicit expressions for the constants 
appearing in these matrices. We use Eq. (|35|l to relate 
the overlined coefficients to the bare coefficients, which 
are defined in Eqs. H33|l and H24|l. We will work to second 
order in the perturbations from isotropic exchange. To 
this order 



Then, keeping only relevant anisotropy corrections, we 
find that the constants in Eq. H43|l and (|44ll are 



<i) 



sin(20) = § = 



and 



cos(26l) 



±±^ I ±±y 



2m ' 



(Cl) 



(C2) 



where Hy and Hz are given in Eq. (|31|l . so that 
sin(2(?) = (—^ " ' 



where 



Then 



3Ji 



Ji 



(C3) 



VSDz Cy 



S{2e) = l~^iJyZ-Cyz-V3Dyf 



(C4) 



(C5) 



Equation P5|l gives 

A = Cxx + C^vv ~ Cxx + ^yy 2Czz "iCyz sin(20) 
= A - ■iCyz sin(26l) . (C6) 



= A + 2E^^^ + 2E^^^ + [2E^^^ - 2E^^^] cos(26') 

-[4i?W+3C,,]sin(20) 
3 1 

3 1 r- 



2[jyz — Cyz ~ ^/^Dy 

9J2 



1 



AViDy - AJyz + 6Cyz 



JyZ Cyz \/3Dy 



2J + A + ^Aj + ^r]j - 2V3Dz 



3Ji 



[jyZ C^yZ '^/^Dy][jyz ^sfiDy 2Cyz] 



-E^^^ ~E 

^x -^y 



{1) 



-i?«-[£;W+i?W]/2 
-^[4'^-i?f^]cos(20)-4'^sin(20) 



V3, 



— ^ J XX ^ Jyy + 2 g '^^ 



T Jyy c^J zz 



V3 



Dz 



J ZZ r> 'J.. 



XX g'Jyy ^ 



2 [J yz Cyz 



~[V3Dy - Jyz][JyZ - Cyz - V3Dy] 



It 3 11^ V3„ 



8 ' 24 



1 

Bq ~ C'xx ~ C 

- C -i 

— ^xx ^ 



Jyz ^ Cyz — V^Dy 



— \/3Dy + Jyz ~\~ Cyz 



Cyy ~\~ C zz 



C zz Cyy 



cos(26') 



+Cyz sin(26') 

Cxx Cyy 



2Cyz 



[Jyz ~ Cyz — V^Dy] , 



-E'x'> + Ey 



^(1) + ^;(i) 

y ' z 



sin(20) 
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7 . 1 



cos(26') 
V3. 



^Jxx gJyy 2'-^^^- ~^ 4 



3 1 



' 2[Jy. - Cy. - V3Dy]^ 

f ^[x/SDy - Jyz][Jyz - Cy^ - ViDy] 

3.5 1 ^ 



Then we will need the frequency and wave vector depen- 
dent susceptibility given by 

x"^{co,ci) = ((5"(q);5'5(-q))). (D2) 

in terms of which the dynamical structure factor is 



where n{u)) = [e^^ — 1] ^ . We have that 
^Qx"''(c^-*0+,q) 



[jyZ — Cyz — VS-Dy] 



6J 



,5(6^ - E„, + £;„)(n|5"(q)|m)(TO|5^(-q)|n) 



= 2^^^ = 24^) cos 61 + 24^) sin 6* 

= 24^) + sin(26') 

= -{Dy + \/iJyz) cos 9 

+ I -v3Ji - Z?. + — I sn 



+ - i;„)(n|5''(-q)|m)(m|5"(q)|n) 



ri(a;) 



But using Eq. (|C3|) this gives 

Jyz Cyz \/3-^y 



3Ji 



- V,,8{lo - + En){n\S"{q)\m){m\S^{-q)\n) 

7n.n \, 

(C7) -pr^Siuj + Er, - Era){m\S^{^q)\n){n\S''iq)\m)'^n{uj) 

^ Y,Pn5{oJ ~ Ern + E„)[l ~ e-^(^'"-^")]n(c^) 

x(n|5"(q)|m)(m|5''(-q)|n) 
1 „ 



Ji 



%/3A,7 V3?7j 
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p/3(£;,„-£;,.) _ 1 



— Jy. + + 



x(n|5"(q)|m)(m|5''(-q)|n) 
^ P™<5(l^ + (n|5" (q) |m) (m|5'' (-q) |n) . (D4) 



Aj - 77J - 2%/3i:'2 + Cy_x + 3 J2 ■ (C8) We write 



Also 



a(2) = 24'^cos0 + 24') sin ^ 
= (-V3J2)sine' 
V3J2 



3Ji 



4--Cy^-v^-Dy ■ (C9) 



APPENDIX D: DYNAMICAL STRUCTURE 
FACTOR 

We first define the Green's function 

(ji\A\ra){m\B\n) 



T,R,P 

where the subscript £ indicates a component with respect 
to the local canted axes. Then, to linear order in the Bose 
operators we have 



W — Em + E. 

n\B\m) {m\A\n 
w + E„i — En 



T,R 



n^,iT)s^{R,T) + n^y{T)sUR,T) 



(Dl) 



27V ^ 



,iq-(T+R) 



T,R 
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Now set 



7^„,(T)[a(R,r)+at(R,x) 



in^yiT)[a(R,T) - a\R,T)] . (D6) 



Fair) = Uax{T) + iUayir) 



(D7) 



E 



rc(r)4^'(q)* -r„(r)*4^'(q)- ] xUq) 



+ r„(r)-4^'(q)-r<,(r)4^"(q) U^(-q) 



For general reciprocal lattice vector one has 



■ (D8) 



Then 
5°(q) 



^ ^>q-(R+-r) |^r„(T)a+(Rr) + T,, (r) * a(RT)] 



IT. 



r„(r)Qt(q)+r„(r)-a^(-q) 



rc(r) ct^'(q)*xt(q)-4^'(q)X^(-q) 



+ro(r)- 4"'(-q)x^(-q) - dt""' (-q)*^^ (q) 



5"(q + G) = ./S/2 



E 



+ r4r)*4^'(q) -r4r)d(.'"(q) U„(-q) 



TD9) 



This result indicates that the scattering intensity will 
oscillate when e**^ ''' changes sign as one goes from one 
Brillouin zone to the next. 
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